LIGHT SCATTERING FROM MESOSCOPIC OBJECTS 
IN DIFFUSIVE MEDIA 



by J.M. Luck^' a ) and Th.M. Nieuwenhuizen( 2 - b ) 

(1) C.E.A. Saclay, Service de Physique Theorique, 91191 Gif-sur-Yvette cedex, France. 

(2) Van der Waals-Zeeman Laboratorium, Valckenierstraat 65, 1018 XE Amsterdam, The 
Netherlands. 



Abstract. No direct imaging is possible in turbid media, where light propagates dif- 
fusively over length scales larger than the mean free path I. The diffuse intensity is, 
however, sensitive to the presence of any kind of object embedded in the medium, e.g. 
obstacles or defects. The long-ranged effects of isolated objects in an otherwise homo- 
geneous, non-absorbing medium can be described by a stationary diffusion equation. In 
analogy with electrostatics, the influence of a single embedded object on the intensity 
field is parametrized in terms of a multipole expansion. An absorbing object is chiefly 
characterized by a negative charge, while the leading effect of a non-absorbing object is 
due to its dipole moment. The associated intrinsic characteristics of the object are its 
capacitance Q or its effective radius R e s, and its polarizability P. These quantities can be 
evaluated within the diffusion approximation for large enough objects. The situation of 
mesoscopic objects, with a size comparable to the mean free path, requires a more careful 
treatment, for which the appropriate framework is provided by radiative transfer theory. 
This formalism is worked out in detail, in the case of spherical and cylindrical objects of 
radius R, of the following kinds: (i) totally absorbing (black), (ii) transparent, (iii) totally 
reflecting. The capacitance, effective radius, and polarizability of these objects differ from 
the predictions of the diffusion approximation by a size factor, which only depends on the 
ratio R/£. The analytic form of the size factors is derived for small and large objects, 
while accurate numerical results are obtained for objects of intermediate size (R ~ £) . For 
cases (i) and (ii) the size factor is smaller than one and monotonically increasing with R/£, 
while for case (iii) it is larger than one and decreasing with R/i. 



P.A.C.S.: 42.30.-d, 42.25.Bs, 42.68.Ay 

(a) e-mail: luck@spht.saclay.cea.fr 

(b) e-mail: nieuwenh@phys.uva.nl 



1 



1 INTRODUCTION 



Light undergoes multiple scattering when propagating through an inhomogeneous 
medium over distances much larger than one mean free path (see ref. [1] for a review). 
Objects embedded in such a turbid medium are invisible from the outside: no ballistic 
imaging is possible. Obstacles or defects can, however, be detected by a careful analysis of 
the perturbations of the diffuse intensity field. Imaging through turbid media has been the 
subject of many investigations. It has a wide range of applications, especially in the field 
of medical imaging. Just to mention one example, at this moment the Philips Laboratories 
in Eindhoven (the Netherlands) are testing prototype devices to detect breast cancer by 
diffuse light propagation. This method is non-invasive: as opposed to conventional X-ray 
tomography, it does not cause radiation damage, even if applied to a large population over 
a long period of time. 

In ref. [2] it has been shown theoretically, and demonstrated experimentally, that the 
diffuse image of an object embedded in a homogeneous turbid medium, whether absorbing 
or not, can be quantitatively described by the diffusion approximation. This framework, to 
be recalled in section 2, is advantageously rephrased in the language of electrostatics. The 
long-distance effect of an embedded object on the diffuse intensity field /(x) is characterized 
by a multipole expansion. If the object absorbs radiation, the leading term is a negative 
charge q, while a non-absorbing object is chiefly described by a dipole moment p. Like in 
electrostatics, higher-order multipoles are most often negligible. This framework leads to a 
quantitative prediction of the diffuse image of the object, i.e., the profile of diffuse intensity 
transmitted through a thick slab near the embedded object. In ref. [2] the charge and the 
dipole moment of a sphere and of a cylinder have been calculated within the diffusion 
approximation, namely by solving the diffusion equation in the presence of a spherical or 
cylindrical object. This approach is justified a priori for macroscopic objects, with a size 
much larger than the mean free path I. Further results on imaging of objects within the 
framework of the diffusion approximation can be found in refs. [3-6]. 

The diffusion approximation breaks down in the case of mesoscopic embedded objects, 
with a size comparable to the mean free path. Recently, Lancaster and Nieuwenhuizen 
have investigated the cases of point scatterers and of small mesoscopic scatterers [7]. In 
the situation of point scatterers, the analysis is rather simple, although short-distance 
singularities have to be regularized, introducing thus new effective parameters. A partial- 
wave expansion was introduced in the more difficult case of extended objects. 

The goal of the present work is to propose a more general approach to the quan- 
titative calculation of the charge and dipole of a mesoscopic object. We shall put the 
emphasis on the associated intrinsic characteristics: capacitance Q, effective radius -R e ff 5 
polarizability P. The natural framework to evaluate these parameters is radiative transfer 
theory (RTT) [8-11], at least in the semi-classical regime, corresponding to most exper- 
imental situations, defined by I ^> A, with A being the wavelength of radiation. RTT 
provides a quantitative description of the intensity over mesoscopic length scales, down to 
the mean free path £, whereas the diffusion approximation only describes the variations of 
the diffuse intensity I(x) over distances much larger than I. This advantage of RTT with 
respect to the diffusion approximation, especially when dealing with skin-layer phenomena 
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near boundaries, has been emphasized in several recent works [12-14]. For definiteness 
and simplicity, we consider multiple isotropic scattering of scalar waves in a non-absorbing 
medium. We focus our attention on spherical or cylindrical embedded objects, with radius 
R, and restrict the analysis to the following three cases: 

(i) totally absorbing (black) object, 

(ii) transparent object, with the same optical index as the medium, 

(iii) totally reflecting object. 

In section 2 we recall the diffusion approach of ref. [2] . In particular we express the charge 
and the dipole of an embedded object in terms of intrinsic characteristics, its capacitance 
Q or its effective radius R e s, and its polar izability P, and we recall the predictions of 
the diffusion approximation in the cases listed above. In section 3 we employ RTT in 
order to derive the capacitance Q or the polarizability P of a sphere in the same three 
cases. Analytical predictions are obtained for small spheres (R <C €) and large spheres 
(R ^> £), while the appropriate RTT equation is solved numerically in the general case 
(R ~ £). The same analysis is then performed in section 4 for the effective radius _R e ff or 
the polarizability P of cylindrical objects. Section 5 contains a brief discussion. 

2 GENERAL FRAMEWORK 

2.1 Multipole expansion 

Our goal is to investigate the diffuse image of a small object embedded in a homoge- 
neous, non-absorbing multiple-scattering medium. The size L of the sample is assumed to 
be very large with respect to the mean free path £, so that the diffusion approximation can 
be used to describe the intensity /(x), everywhere except in a skin layer around the object 
under study, with a thickness of order I. The intensity J(x) is therefore the solution of the 
stationary diffusion equation 

V 2 /(x) = 0, (2.1) 

with appropriate boundary conditions, taking account of sources, of the geometry of the 
sample, and of the presence of the embedded object. Let /o( x ) be the solution of eq. (2.1) 
in the absence of the object. The total intensity reads 

7(x) = J (x) + 5I(x). (2.2) 

The main goal of the present work is to provide a quantitative characterization of the 
disturbance 8I(x) of the diffuse intensity, representing the diffuse image of the embedded 
object. 

• The case of a sphere 

Consider first the case where the embedded object is a small sphere, with radius 
R, located at the position xo in the medium. It is expected on general grounds that 
the disturbance <5/(x) can be characterized by a multipole expansion, in analogy with 
electrostatics [2, 4-7]. It turns out that only the first two of them, namely the charge 
q and the dipole moment p, play a role, while higher-order multipoles are most often 
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negligible. Assume for a while that the medium is infinite. Setting r = x — x and r = |r|, 
the disturbance of the intensity far enough from the object assumes the form 

8I(x) = ( q + p.V) 1 -= q -- 1 ^. (2.3) 

This expression, which is a solution of the diffusion equation (2.1), holds as soon as the 
distance to the object is large with respect to the mean free path i. 

The linearity of the problem and the isotropy of the spherical defect imply that the 
charge q is proportional to the local intensity io( x o) i n the absence of the object, while 
the dipole moment p is proportional to its gradient. We therefore set 

g = -Q/ (xo), p = -PVI (x ), (2.4) 

and we propose to call Q the capacitance of the spherical object, and P its polarizability. 
These two parameters are intrinsic characteristics of the object. Their dependence on its 
radius and on its optical properties will be investigated in detail in the following. 

The capacitance Q is non-zero (and in fact positive) only if the sphere is absorbing. 
Indeed, eq. (2.3) implies that the total flux of radiation absorbed by the object reads 

5F = -Anq = 4tvQI(x ). (2.5) 

For a non-absorbing object, the leading contribution is that of its dipole p. The corre- 
sponding polarizability P can be either positive or negative, depending on the nature of 
the sphere. We finally observe that Q has the dimension of a length, while P has the 
dimension of a volume. As a consequence, for a sphere of radius R, these parameters 
scale as Q ~ R and P ~ i? 3 , roughly speaking. More accurate predictions will be derived 
throughout the rest of this paper. 

• The case of a cylinder 

The present formalism can be extended to other objects. Let us consider the case 
of a rod, modeled as an infinitely long cylinder, with radius R, whose axis is parallel to 
the z-axis. Assume that the sample is infinitely long in this direction. Then both io( x ) 
and the disturbance SI(x.) only depend on the two-dimensional perpendicular component 
xj_ = (x, y). Eq. (2.3) is replaced by 

W(xx) = (q + P • V) In r f = qln ^ - ^ (2.6) 

where = x^ — (xo)j_, so that r = |rj_| is the distance to the axis of the cylinder. 
The length scale ro, requested by dimensional analysis, is determined by conditions at 
the boundaries of the sample, so that it is in general proportional to the sample size. 
This sensitivity of SI(x) to global properties of the sample arises because the logarithmic 
potential of a point charge in two-dimensional electrostatics is divergent at long distances. 
As a consequence, the capacitance Q is not intrinsic to the cylinder. An intrinsic quantity 
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is its effective radius R e s, denned by the condition that the cylindrically symmetric part of 
the total intensity, including the charge term in the right-hand side of eq. (2.6), vanishes at 
a distance r = _R e ff from the axis of the cylinder. We thus have Io( x o) + Q hi( r o/-Reff) = 0, 
hence 

Q = -^nr- (2-7) 

In 

-Reff 

The general remarks of the previous paragraph also apply to a cylinder. We observe that 
i? e ff and ro have the dimension of a length, while the capacitance Q is dimensionless, 
and the polarizability P has the dimension of an area. As a consequence, the orders of 
magnitude i? e ff ~ R and P ~ R 2 can be expected for a cylinder of radius R. More accurate 
predictions will be derived throughout the rest of this paper. 

2.2 Transmission in a thick-slab geometry 

In order to illustrate the definitions of the capacitance Q, of the effective radius R e s, 
and of the polarizability P given in section 2.1, we recall the derivation of the diffuse image 
of a small object, i.e., the profile of transmitted intensity, through a thick-slab sample. This 
geometry is often encountered experimentally [2]. 

We consider an optically thick slab (0 < x < L), with L ^> £, assumed to be infinite 
in the two transverse directions. The left side of the sample (x = 0) is subjected to an 
incident plane wave, so that the solution of eq. (2.1) when the object is absent reads 

Jo(x) = Jo (l - \ ) • (2.8) 

The intensity T(y, z) transmitted through the sample, and emitted on the right side 
(x = L) at the point (y, z), is proportional to the normal derivative of the diffuse intensity 
at that point: 

T(y,z) = -Ki dI{X ^ Z) 



9X x=L 

When there is no embedded object, eq. (2.8) yields a uniform transmission 



(2.9) 



T = ™ (2.10) 

The prefactor K bears the dependence of the transmitted intensity on the angle and 
polarization of the incident and emitted radiation. If these characteristics are not resolved, 
namely if the "all-in, all-out" intensity is measured, the prefactor K assumes the value 
K = 47r/(3A 2 ), with A being the wavelength of radiation. This universal result is analogous 
to the Boltzmann formula for the conductivity [12, 1]. The dependence of K on incidence 
angles and polarizations is non-trivial, as it involves microscopic characteristics of the 
medium, such as the anisotropy of the scatterers. These features cannot be accurately 
described within the diffusion approximation. Indeed, they originate in the presence of 
skin layers near the boundaries of the sample, with a thickness comparable to i, where 
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the free radiation is converted into diffuse intensity, and vice versa. The appropriate tool 
to investigate skin-layer phenomena is the Schwarzschild-Milne equation of RTT [12-14]. 
Ref. [1] contains an overview of these recent works. 

• The case of a sphere 

Consider first the case where the embedded object is a sphere, with capacitance Q 
and polarizability P, located at x = xq, y = z = 0. Its charge and dipole moment read 

Pla 



with e x being the unit vector of the positive x-axis. 

Along the lines of ref. [2] , the total intensity can be calculated by summing the infinite- 
space expression (2.3) over a periodic double array of images. Setting p = (y 2 + z 2 ) 1 ^ 2 , we 
obtain 



+oo 



/(x) = J (l-|)+g £ 



1 



[(a; - x + 2nLY + p 2 ] 1/2 [(a; + x + 2nL^ + p 2 ] h 2 



n= — oo 
+oo 



Ei x — xq + 2nL x + xq + 2nL 



n= — oo 



[(a; -x + 2nLY + p 2 ] 3/2 [(x + x + 2nLf + p 2 ] 3/2 

(2.12) 

The transmitted intensity T(y, z) = T(p) then reads 



T(p) = T [l-2Q(L-xo) £ L ""° + 2nL 

n _oo [(L -x + 2nL) 2 + p 2 ] 

-IP 2 " X ° + 2UL)2 ~ ^ 



3/2 



(2.13) 



[(L -x + 2nL) 2 + p 2 ] 



5/2 



The diffuse image of the embedded sphere thus extends over a region of transversal 
size p ~ L — xq. If the sphere is absorbing, the leading contribution of its charge is negative, 
with a relative magnitude proportional to Q/L, i.e., to R/L. If the sphere is not absorbing, 
the leading contribution of its dipole is proportional to P/L 3 , i.e., to R 3 /L 3 , in magnitude. 
It is negative for P > (e.g. for a reflecting sphere), and positive for P < (e.g. for a 
transparent sphere). 

• The case of a cylinder 

Consider the case where the embedded object is a cylinder, with given effective radius 
i? e ff and polarizability P, whose axis is parallel to the z-axis, and contains the point x = xq, 
y = 0. The total intensity can be calculated by applying the method of images to eq. (2.6). 
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We thus obtain 



-|-oo 



(x — x + 2nL) 2 + y 2 
(x + xn + 2nL) 2 + v 2 

n= — oo 



E( x — x + 2nL x + xq + 2nL 



(x - x + 2nL) 2 + y 2 (x + x + 2nL) 2 + y 2 
This result has the following alternative closed-form expression 

n(C + x ) 



(2.14) 



Ti \ T ( ~\ X \ T> (\ ■ 7r(C-^0) , . 

I (x) = I 1 1 — — ) — q Re I m sin — — - — - — In sin 



2L 



n 



n((-x ) n(( + x ) 



~ P 2L Re [ COt ^I— + COt 2L 



(2.15) 



in terms of the complex variable £ = x + iy. 

The charge q can be calculated by expressing that the first line of the right-hand 
side of eq. (2.15), including the free term and the charge term, vanishes at the distance 
R e f£ <C L from the axis of the cylinder. We thus obtain 



q = —QIo ( 1 - -r) > witn Q = rvT" an d r o = — sin (2.16) 

\ L / l n U 7T L 

-Reff 



in agreement with the general results (2.4), (2.7). This example shows that the length 
scale ro is proportional to the sample thickness L, and that it depends on the position of 
the embedded object. 

Finally, the transmitted intensity T(y, z) = T(y) reads 



T(y) = T 



( . Kx , , 7tx Try \ 

■ « n V LJ ^'El.^^V L 2 ( « , Try\ 2 



\ v ^cos^ + cosh^ ^ ( cos t« + cosh my j 

(2.17) 

The diffuse image of the cylinder again extends over a region of transversal size y ~ 
L — xq. If the cylinder is absorbing, the leading contribution of its charge is negative. Its 
relative magnitude, proportional to Q, depends logarithmically on the sample thickness L 
via ro, according to eq. (2.16). If the cylinder is not absorbing, the leading contribution of 
its dipole has a magnitude of order P/L 2 , i.e., R 2 /L 2 . 

The characteristic profiles of transmitted radiation recalled above have been observed 
experimentally by den Outer et al. [2], in the case of pencils (black cylinders) and glass 
fibers (transparent cylinders with a non-trivial index mismatch). 
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2.3 Diffusion approximation 

When the radius R of the embedded object is much larger than £, its characteristics 
Q or i? e ff and P can be calculated within the diffusion approximation, along the lines of 
ref. [2]. Indeed, the diffusion approximation only breaks down in a skin layer of thickness 
I < i? around the surface of the object. We summarize, for completeness, the main 
predictions of this approach, in the three cases listed in the Introduction. 

2.3.1 The capacitance of a totally absorbing (black) object 

• The case of a sphere 

Consider first a totally absorbing sphere of radius R ^> £, located at the origin in 
an infinite medium. Within the diffusion approximation, the intensity i"(x) = I(r) is a 
solution of eq. (2.1), more precisely a spherically symmetric harmonic function, with the 
absorbing boundary condition I(R) = at the surface of the sphere. A basis of harmonic 
functions is given by {1, 1/r}, so that the required solution is proportional to 

Ur) = 1 - -. (2.18) 

r 

By identifying this expression with eqs. (2.3), (2.4), we get the simple result 

Q diff = R (2.19) 



for the capacitance of an absorbing sphere, in the diffusion approximation. We recall that 
this result is expected to hold in the regime R ^> t. The dependence of Q on the ratio Rjt 
will be investigated in section 3.2. 

• The case of a cylinder 

Consider a totally absorbing cylinder of radius R ^> £, whose axis is parallel to the z- 
axis and passes through the origin. A basis of cylindrically symmetric harmonic functions 
is given by {1, lnr}, so that the required solution is proportional to 

I(r)= In j. (2.20) 

We thus have 

R d J = R- (2.21) 

The effective radius of an absorbing cylinder coincides with its geometrical radius, within 
the diffusion approximation. Again, this result is expected to hold for R 3> £ only, while 
the dependence of _R e ff on Rjt will be investigated in section 4.2. 
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2.3.2 The polarizability of a non-absorbing object 

• The case of a sphere 

Consider a non-absorbing sphere of radius R ^> £, located at the origin in an infinite 
medium. Within the diffusion approximation, the only property of a material which is 
naturally involved is the diffusion constant of radiation in its bulk. So, assume that the 
diffusion constant reads D\ in the medium and D 2 in the sphere. 

Let the intensity be ii(x) outside the sphere (r > R) and I 2 (x) in the sphere (r < R). 
These functions obey eq. (2.1), together with the conditions 

h = h, D 1 n-VI 1 =D 2 n-VI 2 (r = R), (2.22) 

expressing respectively the continuity of the intensity and the conservation of the flux at 
the boundary of the sphere. Here n denotes the outward normal to the sphere. 

We look for a solution to eq. (2.1) which behaves as i"(x) ~ x far from the sphere. 
This formally amounts to calculating the dipole moment of a dielectric sphere in a uniform 
electric field. It can be checked that the solution assumes the form 




J 2 (x) = Ax. (2.23) 



The outer solution ii(x) is the superposition of a uniform gradient and of a dipole field, 
while the inner solution h(x.) is a uniform gradient. The dipole component of ii(x) yields 
the polarizability P, again by eqs. (2.3), (2.4). 

The conditions (2.22) yield two linear equations for the unknowns P and A, hence 



2D 1 + D 2 



R . (2.24) 



In the present work we are mostly interested in the limit D 2 — > 00, which yields the 
polarizability of a transparent sphere, with the same optical index as the medium, 

pdiff = _ jR 3 ? (2.25) 
and in the limit D 2 — * 0, which yields the polarizability of a totally reflecting sphere, 

pdiff = ^_ ( 22 g) 

Again, these results are expected to hold in the regime R ^> I. The dependence of both 
polarizabilities on Rjl will be investigated in sections 3.3 and 3.4. 
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• The case of a cylinder 

Consider a non-absorbing cylinder of radius R 3> £, located as above. Assume again 
that the diffusion constant reads Di in the medium and D 2 in the cylinder. The intensity 
field can be calculated along the lines of the above case of a sphere. We thus obtain 

Di - Do o , 

P = — -R 2 . 2.27 

D 1 + D 2 V ' 

Again, the limit D 2 — > oo yields the polarizability of a transparent cylinder, with the same 
optical index as the medium, 

pdiff = _ R 2^ 

while the limit D 2 — > yields the polarizability of a reflecting cylinder, 

pdiff = R 2 ( 229 ) 

The dependence of both polarizabilities on R/£ will be investigated in sections 4.3 and 4.4. 

3 RADIATIVE TRANSFER THEORY: THE CASE OF A SPHERE 

We turn to the calculation of the capacitance Q or the polarizability P of mesoscopic 
spheres, with a radius R comparable to the mean free path £, in the three cases listed in the 
Introduction. As recalled there, the appropriate framework is RTT, at least in the regime 
£ ^> A, corresponding to most experiments. In the following, we consider for simplicity 
isotropic scattering of scalar waves. 

Throughout sections 3 and 4, all lengths are measured in units of the mean free path 
£, for convenience, unless otherwise stated. So R stands for the dimensionless ratio R/£, 
and so on. 

3.1 Basics of radiative transfer theory 

We first recall the basic concepts of RTT, following the book by Chandrasekhar [8]. 
In a time-independent situation, the basic quantity is the specific intensity 7(x, n), namely 
the amount of intensity located at the point x and propagating into the direction of the 
unit vector n. In the case of isotropic scattering, the local form of the RTT equation reads 

n- VI(x,n) = J(x) -/(x,n), (3.1) 

in dimensionless units, where the source function X(x) reads 

!(*) = / *£>/(*,»), (3.2) 

with du;(n) being the element of solid angle around the direction n. The source function 
J(x) will be identified with the diffuse intensity J(x) involved in the diffusion approxima- 
tion. 
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The RTT equation (3.1) admits the formal solution 



/•OO 

J(x,n)= / e- s dsJ(y(x,n, S )), (3.3) 
Jo 

where y(x, n, s) is the point such that there is a ray of optical length s, starting from point 
y and arriving at point x while pointing into direction n. We thus obtain the following 
linear integral form of the RTT equation for the source function 

X(x) = J ^l^°°e- s d S J(y(x,n, S )), (3.4) 

which will be the starting point of all subsequent developments. 
In an infinite medium, without any object, we have 

y(x, n, s) = x - ns, (3.5) 

so that eq. (3.4) can be recast as a convolution integral equation of the form 

J(x) - J d 3 yM(x - y)J(y) = 0, (3.6) 

with a convolution kernel 

M(*- y ) = /^V^>( X - y -n s ) = i ^. (3.7) 

Its Fourier transform reads 

M (q) = / M(x) e ^x = / = = (3.8) 

with q = |q|. Therefore the left-hand side of eq. (3.6) reads approximately (— 1/3) V 2 X(x), 
up to higher-order derivatives. This demonstrates that RTT incorporates the diffusion 
approximation in the regime where distances are large with respect to the mean free path 
£, corresponding to q <C 1 in reduced units. 

3.2 The capacitance of a totally absorbing (black) sphere 

By definition of a totally absorbing object, any ray which hits it corresponds to radi- 
ation which is absorbed. A typical such ray, to be discarded from the integral in eq. (3.4), 
is shown as TX in Figure 1. 

We introduce spherical co-ordinates (r, 9, ip) , such that (x,y,z) = (r cos 9, r sin 9 cos tp, 
r sin 9 sin ip) . Because of the rotational invariance of the problem, we can choose to locate 
the observation point X on the positive x-axis, at a distance r > R, and the ray in the 
x-y plane, which coincides with the plane of Figure 1. The direction n is marked by the 
angle a (0 < a < 7r), while the incidence angle of the ray at the surface of the sphere is f3 
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(0 < (3 < 7r/2). Some useful distances and other quantities related to Figure 1 are listed 
in Table 1. 

We are thus led to split the range of integration over n and s in eq. (3.4), and in all 
similar integrals to appear in the following, into two domains, respectively called / and 
II, corresponding to rays which hit the sphere and rays which do not. Table 1 also gives 
analytical definitions of these domains, and of the following shorthand notations 

for the integrals over either domain, and for the full integral over the whole range of 
parameters (0 < a < rc and < s < oo). The latter is normalized: 

1 = 1. (3.10) 

full 

In the following, it will be advantageous to change variable from a to (3 in integrals over 
the domain /, and to use the second form of the integration measure given in Table 1. 

The RTT equation (3.4) for a totally absorbing sphere thus reads 



J(x) = // J(x-ns). (3.11) 



II 



Along the lines of section 2.3.1, we look for a rotationally invariant solution J(x) = J(r) 
to eq. (3.11), going to unity at large distances. Eq. (3.11) for the function J(r) reads 



J(r) = JJ^Jip), (3.12) 



with 

p= |x-ns| = (r 2 + s 2 - 2rscosa) 1/2 . (3.13) 

We set 

J(r) = l-j + F(r), (3.14) 

where F(r) represents the short-ranged correction to the solution (2.3) of the diffusion 
equation. This function therefore describes the skin layer around the sphere. It is expected 
to fall off exponentially, with a decay length given by the mean free path, i.e., F(r) ~ e~ r 
in reduced units. Some algebra allows to recast eq. (3.12) as an equation for F(r) itself: 

F(r) - [[ F{p) + [[ F{p) = S 1 (r)+QS 2 (r), (3.15) 

''full JJl 

with 

S 1 {r) = -J]l, S 2 (r)=f(r)+JJ i - p , (3.16) 
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and 

r i/full P r Jr S 

• Small-radius behavior 

The behavior of the capacitance Q for < 1 can be derived by taking the Fourier 
transform at q — > of eq. (3.15), or equivalently its integral over all space. Indeed, 
the first two terms of the left-hand side cancel out by this procedure, because of the 
behavior (3.8) of the Fourier-transformed convolution kernel M(q). The space integral 
of f(r) is elementary, while the other contribution to ^(r) is negligible. The integral of 
S± (r) can be performed exactly, by changing variables in domain I from a to and from 
R < r < oo to Rcos/3 < w < oo, according to the definitions given in Table 1. We thus 
obtain 



"OO -1 

,2 - -' ^ 1 



L 



f 

/ r 2 dr/(r) = -, 

oo t->2 /■7r/2 />oo p2 

(3.18) 



t->2 /*vr/ z /.oo nk 

r 2 dr^i(r) = -— / sin /3 cos /3 d/3 / e flco8 ^- w d«; = -=j 

2 Jo J Rcos/3 4 



i? ^ </0 J Rcos/3 4 

By inserting these expressions into eq. (3.15), we obtain the following prediction 

O E>2 

Q*^. (3-19) 

in physical units, for R <ti I. The leading correction to this result is of relative order R/£. 
It cannot be predicted by elementary means, since it involves, among other contributions, 
the space integral of the third term in the left-hand side of eq. (3.15). 

The power of Rjt involved in the leading correction to the estimate (3.19) can be 
found more rapidly as follows. Consider the integral (3.18) of the source function S\{f). 
The lower bound of this integral has been set to the natural value r = R, since domains 
I and II are well-defined for r > R only. Changing the lower bound from r = R to 
r = 0, say, would alter the value of the integral by an amount of relative order R/£. By 
consistency of the whole approach, the relative correction to eq. (3.19) is expected to be 
of the same order. This rule of thumb will be used throughout the following. 

The leading-order result (3.19) has been obtained independently in ref. [7]. Paasschens 
and 't Hooft [6] obtain the estimate Q ~ R 2 /£, by evaluating the average dwell time, i.e., 
the mean time spent by radiation inside the object. The discrepancy between their result 
and eq. (3.19) can be explained by noticing that the appropriate quantity to consider in 
their approach is rather the inverse of the mean inverse dwell time [15]. 

• Large-radius behavior 

If the radius of the sphere is large (R 3> 1), its surface is flat at the scale of the mean 
free path, where F(r) exhibits appreciable variations. Hence in this limit eq. (3.12) is 
expected to become a RTT equation in a half-space. This is indeed the case. Let us set 

r = R + T (0<r<oo), (3.20) 
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and J(r) = T(r). We have, to leading order as R — > oo, 

p = R + r', with t' = t — s cos a, (3-21) 
this quantity being positive in domain II, so that eq. (3.12) can be recast as 

/•oo 

r(r)- / dr / M 1 (r-r / )r(r / ) =0, (3.22) 



where 



M 1 (t)= If 5(r-scosa)= [ ^ e -M/^ = -i E i( - |r|), (3.23) 

i/full JO 2 ^ 2 

with Ei(x) being the exponential integral function. The full integral appears in eq. (3.23), 
while domain / does not contribute, because cto = tv/2 in this regime. 

Eq. (3.22) is nothing but the homogeneous Schwarzschild-Milne (SM) integral equa- 
tion of RTT in a half-space with a free boundary. The function M\(x) is the associated 
Milne kernel for isotropic scattering, with the subscript 1 reminding that M\(x) is a one- 
dimensional projection of the RTT kernel M(x) of eq. (3.7). The normalized solution to 
eq. (3.22) behaves as [12-14, 1] 

T(t)^t + t (r»l), (3.24) 

where 

1 r n / 2 Ht tan 2 t 

r = - / In tan x = 0.710446090 (3.25) 

7T Jo sin x 3(1 - x cot x) 

is Milne's extrapolation length [8, 11]. This term reminds that the extrapolation of the 
diffusive behavior (3.24) vanishes at r = — r . 

By requiring that expression (3.24) is proportional to the solution of the diffusion 
equation, J(r) = 1 — Q/r, in the range 1 C r < J? [see eq. (3.14)], we obtain the estimate 

Q ~ -FT—-0 ~ # f 1 " ^TrV ( 3 -26) 



in physical units, for R ^> £. The leading-order behavior of eq. (3.26) agrees with the 
prediction (2.19) of the diffusion approximation, Q dlff = R. The first correction to this 
result has a simple interpretation: the diffusion approximation has to be corrected by 
taking into account a skin layer of thickness tqI, and requiring that the diffusive field does 
not vanish exactly at the surface of the sphere, but at a distance Toi inside the sphere. 

For arbitrary values of the radius R of the sphere, its capacitance Q is such that the 
inhomogeneous linear equation (3.15) for F(r) has a solution which decays rapidly as r — > 
oo. This algorithm can be implemented numerically, along the lines of section V of ref. [12]. 
The idea is to discretize eq. (3.15) into an inhomogeneous system of linear equations, and 
to solve the latter system numerically, first with a right-hand side equal to Si(r), yielding 
a solution -Fi(r), then with a right-hand side equal to S2 (r), yielding a solution -F 2 (r), 
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so that Q is determined by requiring that the linear combination F(r) = Xi(r) + QF 2 (r) 
decays rapidly as r — > oo. 

Figure 2 shows a plot of the size factor Q/Q dlf£ against R/t. The numerical data ob- 
tained by the above method smoothly interpolate between the asymptotic behaviors (3.19) 
and (3.26), listed in Table 2. 

3.3 The polarizability of a transparent sphere 

Consider now a transparent sphere, having the same index of refraction as the medium, 
e.g. a cavity containing no scatterers. The status of the rays in domain II, which do not 
hit the sphere, is unchanged with respect to the previous case of an absorbing sphere: 
y(x, n, s) is again given by eq. (3.5). A typical ray in domain I, which hits the sphere, is 
the ray TX of Figure 1. The part BC does not contribute to the optical length of the ray, 
which reads s = \TB\ + \CX\ = \TB\ + sq. The point T is therefore located at 

y T (x,n, s) =x-s T n, with s T - s = \BC\ = 2Rcos/3 = 2(R 2 - r 2 sin 2 a) 1/2 . (3.27) 

The RTT equation (3.4) for a transparent sphere thus reads 

J(x) = JJl{* ~ ™ T ) + JJ X(x - ns). (3.28) 

Eq. (3.28) admits the solution X(x) = 1, showing that a transparent sphere has a van- 
ishing capacitance, as anticipated since it does not absorb radiation. In order to calculate 
its polarizability P, in analogy with the outer solution ii(x) of eq. (2.23), we look for a 
solution to eq. (3.28) of the form 

J(x) = xK(r) = rK(r) cos 9, (3.29) 

with K (r) going to unity at large distances. Eq. (3.28) then becomes the following equation 
for the unknown function K(r): 

rK{r) = JJ (r — st cos a) K(pt) + JJ (r — s cos a)K(p), (3.30) 

where p has been defined in eq. (3.13), and with 

p T = ( r 2 + s 2 T - 2rs T cosa) 1/2 . (3.31) 

It can be checked that eq. (3.30) is indeed independent of the particular choice we made, 
namely to put the observation point on the positive x-axis. 

We now set 

K(r) = l + ^ + G(r), (3.32) 
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where G(r) again represents the short-ranged skin-layer correction to the solution (2.3) of 
the diffusion equation. The equation for G(r) reads 

rG(r) + (s cos a— r)G(p)+ // [(st cos a— r)G(pT) + (r— s cos a)G(p)~\ = Si(r)+PS2(r) : 

i/full JJi 

(3.33) 

with 



Si(r) = // (s — st) cos a = — 2P // cos a cos/?, 

JJ '-r n , ■"' ^ (3-34) 

„ , s e 11 I f — st cos a s cos a — r x 

* (r) = -F + i(— J?— + — 

• Small-radius behavior 

The behavior of the polarizability P for P 1 can be derived by taking the Fourier 
transform at q — > of eq. (3.33), after multiplication by cos#. The Fourier transform of 
the function J(x) defined in eq. (3.29) reads 

= r~ / r drP(r)(singr — gr cosgr) / r drP(r) (g — > 0). (3.35) 

Q Jo 3 Jo 

Therefore both sides of eq. (3.33) have to be integrated with a weight r 3 dr. The integral 
of Si(r) can again be estimated by changing variables from a to j3 and from r to w: 

/2 /.oo /i 

> ■ . / e Rcos ^ w wdw = -R 3 - + 4). 



pOO pTT '/ A p 

/ r 3 drS 1 (r) = -R 3 sin /3 cos 2 /3 d/3 / 
i.R io i.fi 



i? cos (3 



3 4 

(3.36) 

The integral of the first contribution to Si (r) is elementary, while the second contribution 
is negligible. 

We thus obtain the estimate 

+ (3.37) 

in physical units. The rule of thumb exposed in section 3.2 predicts that the leading 
correction to the estimate (3.37) is of relative order R 2 /£ 2 . 

• Large-radius behavior 

When the radius of the sphere is large (P ^> 1), eq. (3.30) again becomes a RTT 
equation in a half-space. With the definition (3.20), we have 



p = R + t' , with t' = t — s cos a > 0, in domain PZ", 
Pt = R + t", with t" = s cos a — r > 0, in domain I. 



(3.38) 



Setting K(r) = T(r), eq. (3.30) can be recast as 

p oo 

r(r)- / dr , [M 1 (r-r , ) + M 1 i (r + r / )]r(r / ) = 0, (3.39) 
Jo 
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where 

M[(x)= ff (1 - 2 cos 2 a) 5(x-s cos a) = f ^(1 - 2 / u 2 ) e -l x l^. (3.40) 
i/fuii Jo 2// 

We thus obtain the SM integral equation of RTT in a half-space with partial reflection at 
the boundary, characterized by an effective intensity reflection coefficient 

R{n) = 1 - 2/j 2 . (3.41) 

The reflection coefficient is negative for a small enough, namely for \i = cos a > l/\/2, i.e., 
a < 7r/4. This observation is less surprising if we remember that we are only investigating 
the contribution to the intensity field of the dipole of the spherical object. This contribution 
is not a positive function, because of its angular dependence, given by eq. (3.29). 

The normalized solution to eq. (3.39) behaves as 

r(r)^r + r i (r > 1). (3.42) 

The extrapolation length roi cannot be calculated analytically. We have evaluated it 
numerically, again along the lines of section V of ref. [12], obtaining thus 

r i ~ 0.4675. (3.43) 

By requiring that expression (3.42) is proportional to the solution the diffusion equa- 
tion, K(r) = 1 + P/r 3 , in the range 1 C r < J? [see eq. (3.32)], we obtain the estimate 

R 1 .._ Ig ,( 1 _^f) < (3.44) 



R + 3t i£ V R 

in physical units, for R ^> I. This result agrees, to leading order, with the prediction (2.25) 
of the diffusion approximation, _P dlff = — R 3 . 

Figure 3 shows a plot of the size factor _p/_p dlff against R/£. The numerical data 
interpolate between the asymptotic behaviors (3.37) and (3.44), listed in Table 2. 

3.4 The polarizability of a reflecting sphere 

Consider now a reflecting sphere. The status of the rays in domain II, which do 
not hit the sphere, is unchanged with respect to the two previous cases. A typical ray in 
domain /, which hits the sphere, is the ray RCX on Figure 1. We have s — sq = \RC\, 
while s = \XC\. The point R is thus located at 

y fl (x, n, a) : ( Xr = r C ° S ? + ( 8 ~ f ^^" a) ' (3.45) 
JHy ' ' ' \ y R = s sina + (s - s ) sm(2/3 - a). v ; 

The RTT equation (3.4) for a reflecting sphere thus reads 

J(x) = JJ l(y R ) + JJ J(x - ns) . (3.46) 
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We again look for a solution to eq. (3.46) of the form 

J(x) = xK(r) = rK(r) cos (3.47) 
with K(r) going to unity at large distances. Eq. (3.46) becomes 

rK(r) = JJ x R K{p R )+ JJ {r-scosa)K{p), (3.48) 

where p has been defined in eq. (3.13), and with 

PR=\y R \ = (x 2 R + y 2 R ) 1 / 2 . (3.49) 

Setting again 

K(r) = l + ^ + G(r), (3.50) 



the equation for G(r) reads 



rG(r)+ // (s cos a-r)G(p)+ // [(r-s cos a)G(p) -x R G(p R )] = S 1 (r)+PS 2 (r), (3.51) 

I/full JJl 

with 

S±(r) = // (scosa — r + xr) =2 // (s — sq) cos,9cos(q; — /?), 

„ . . e I xr s cos a — r N 



Small-radius behavior 



r2 JJi \Pr P 3 



The behavior of the polarizability P for R <C 1 can again be derived along the lines 
of the previous section. The integral 

/•oo / >7r /2 poo 

I r 3 dr Si(r) = R 2 sin /3 cos 2 /3 d/3 / e ficos/3_u '(tw cos /3 + Rsin 2 /3)dw 

J R JO J RcosP ^ gg^ 



leads to the estimate 

R 2 l 

P«^, (3.54) 

in physical units, up to a correction of relative order R/i. This result agrees with the 
leading-order prediction of Lancaster and Nieuwenhuizen [7], P ~ cr tot £/(47r), with the 
total cross-section cr tot assuming its geometrical value a to t = ttR 2 . 
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• Large-radius behavior 

When the radius of the sphere is large, eq. (3.48) again becomes a RTT equation in a 
half-space. Indeed, eqs. (3.38) still hold, to leading order. Setting K{r) = T(r), eq. (3.48) 
can be recast as 

/•OO 

r(r)- / cIt'[Mi(t-t') +M 1 (t + t')]T(t') = 0. (3.55) 
Jo 

We thus obtain the SM integral equation of RTT in a half-space with total reflection at 
the boundary. This equation admits a constant solution, which we normalize by setting 

r(r) = i. 

In order to unravel skin-layer phenomena in the present careful analysis of 

corrections to the RTT equation (3.55) is needed. Let us anticipate that we have to work 
up to order 1/R 3 included, and set 

rM i _l ri ( r ) I r 2(r) T 3 (r) 

r(r) = l + ^ r + ^ r + ^ + ---, (3.56) 

with the normalization condition T(0) = 1, hence Ti(0) = r 2 (0) = 1^3(0) = 0. 

The successive correction terms of eq. (3.56) obey inhomogeneous SM integral equa- 
tions, which can be obtained recursively, by systematically expanding the various quantities 
as power series in 1/R, keeping r and s fixed. For example we have 

r sin 2 a r 2 sin 4 a r 3 sin a(5 — 4 cos 2 a) r 4 

S = 3 D+o 5~ "52 H a 7 H ' t 3 ' 57 ) 

cos ol 2 cos ol H, 2cos°aR z 8 cos 'a R* 

and similar series expansions for (3, xr, pr, and so on. Some algebra yields successively 

/•OO 

]?! (r) - / dr' [M t (r - r') + M t (r + r')] {Ti (r) + r' - r) = 0, (3.58) 
Jo 

whose solution, normalized to Ti(0) = 0, is 

r 1 (r) = -r, (3.59) 

and 

f°° 4 

(r)- / dr / [M 1 (r-r / )+M 1 (r + r / )]r 2 (r / ) = --, (3.60) 
Jo <J 



r 2 

whose normalized solution is 



T 2 (r) = 2r 2 . (3.61) 



The determination of the third correction term is less easy. Anticipating the behavior 
^(t) = k\T + k 2 T 2 + k 3 T 3 + 73(1"), we obtain by consistency ks = —10/3, /c 2 = 0, and 

7s(t)- / dr'[M 1 (r-r') + M 1 (r + r')] l3 (r') 

Jo (3.62) 

= ((fci + 8)r 2 - r^M 1 (r) + ^(k 1 + 2-(k 1 + 6)r - r 2 + r 3 )e" r , 

19 



with 73(0) = 0. A necessary condition for eq. (3.62) to admit a bounded solution is that 
the integral of its right-hand side over the range < r < 00 vanishes. This determines 
k\ = —4/5, so that we have 

r 3 (T)=-yT 3 -^T + 7 3(T). (3.63) 

The calculation of the bounded part 73 (r) cannot be performed analytically. 

By requiring that expression (3.56) is proportional to the solution of the diffusion 
equation, K(r) = 1 + P/r 3 , in the range 1 « r <C i? [see eq. (3.50)], we obtain the 
estimate 

1 + (3-64) 

in physical units, for R ^> £. The leading-order behavior of eq. (3.64) agrees with the 
prediction of the diffusion approximation (2.26), _P dlff = _R 3 /2. The correction term in 
£ 2 /R 2 entirely comes from the term in r in eq. (3.63), while eqs. (3.59), (3.61), and the 
leading term in r 3 in eq. (3.63) just reproduce the expansion of the diffusive solution for 
r < R. 

Figure 4 shows a plot of the size factor _p/_p dlff against R/£. The numerical data 
interpolate between the asymptotic behaviors (3.54) and (3.64), listed in Table 2. 

To close up, we mention that the result (3.64) can be recovered by the following more 
direct, albeit less systematic, alternative route. Assuming the behavior 

p3 

P « — + AR, (3.65) 

and systematically using expansions such as eq. (3.57) for so, we can recast the right-hand 
side of eq. (3.51) as 



1 rl 2 sin a da C°° 
Sl (r)+PS 2 (r)^-2 —z— e- s ds$(T,a,s 

n JO L Jr/cosa 



(3.66) 



with 



$(t, a, s) = 4Acos 2 a + (7 cos 2 a — 3)r 2 + (9 — 17 cos 2 a)rs cos a + (10 cos 2 a — 6)s 2 cos 2 a. 

(3.67) 

Setting G(r) = T(t)/R s , eq. (3.51) implies that T(r) obeys an inhomogeneous SM equa- 
tion, with a right-hand side given by eq. (3.66), without the prefactor 1/R 2 . In order for 
this SM equation to have a bounded solution T(r), it is necessary that the r-integral of 
its right-hand side vanishes. We thus obtain after some algebra A = 3/5, in agreement 
with eq. (3.64). The corresponding solution, normalized by the condition T(0) = 0, reads 
f (r) = (3/2) 73 (r). 
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4 RADIATIVE TRANSFER THEORY: THE CASE OF A CYLINDER 



We turn to the calculation of the effective radius _R e ff or the polarizability P of cylin- 
ders. This situation is of special interest, since it corresponds to the experiments by den 
Outer et al. [2]. We recall that throughout this section lengths are measured in units of 
the mean free path £, unless otherwise stated. 

4.1 Radiative transfer theory in cylindrical geometry 

We are led to consider RTT in the presence of a cylindrical object of radius R, whose 
axis is parallel to the z-axis and passes through the origin. In such a circumstance, the 
source function X(x) is independent of the co-ordinate z. We denote by x = xi = (x, y) 
the position in the plane perpendicular to the z-axis, and we use cylindrical co-ordinates 
(r, 9, z), such that (x,y) = (r cos 6>, r sin 9). The x-y plane again coincides with the plane 
of Figure 1. 

The key ingredient of this analysis will be the two-dimensional projection M 2 (a) of 
the RTT kernel M(x) of eq. (3.7), defined in analogy with the one-dimensional projection 
Mi(t) of eq. (3.23), namely 

M 2 (a) = J e- s dsd(a ~ sn\\) = £^ d(/)e- a/cos ^, (4.1) 

where <p is the angle between the unit vector n and the x-y plane, with < <p < tt/2, 
so that n || = cos<p is the length of the projection of n onto this plane. The last integral 
expression of eq. (4.1) will allow us to calculate the kernel function M 2 (a) numerically, 
and to perform the required analytical calculations involving this function. We mention 
for completeness that its derivative is equal to dM 2 (a) / da = —Ko(a), with Kq being the 
modified Bessel function. 

4.2 The effective radius of a totally absorbing (black) cylinder 

We first investigate the effective radius _R e ff of an absorbing cylinder of radius R. In 
analogy with the spherical geometry dealt with in section 3, we split the range of integration 
over a and a into two domains, / and II, respectively corresponding to rays which hit 
the cylinder and which do not. The definitions of these domains, and of the corresponding 
integration measures, are given in Table 1. 

In analogy with eq. (3.11), the RTT equation for an absorbing cylinder reads 

J(x) = ^J(x-na), (4.2) 

with n being the unit vector of the plane with a direction marked by the angle a. Along the 
lines of sections 2.3.1 and 3.2, we look for a cylindrically symmetric solution J(x) = J(r) 
to eq. (4.2). This equation becomes 

j(r)= JL j(p) ' (4 - 3) 
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with 

p = |x - ncr| = (r 2 + a 2 - 2ra cos a) 1/2 . (4.4) 
We normalize the solution J(r) by setting 

J(r) = In ^- + F(r), (4.5) 

where the function F(r), describing the skin-layer effect, is again expected to decay expo- 
nentially as r — > oo. We recast eq. (4.3) as an equation for F(r) itself: 

F{r) - II F(p) + II F(p) = S 1 (r) + S 2 (r) lni? eff , (4.6) 

i/full JJl 



with 



and 



S 1 (r) = f(r)-£\np, S 2 (r) = jj \ (4.7) 

f(r)= II ln^ = / M 2 (a)da\n-. (4.8) 
i/fuii r J r r 



Small-radius behavior 



The behavior of the effective radius i? e ff for R <C 1 can again be derived by taking 
the Fourier transform at q — > of eq. (4.6), or equivalently its integral over all space. The 
space integrals of f(r) and S 2 (r) can be evaluated by changing variables in domain / from 
a to (3, and from R < r < oo to Rcos(3 < w < oo, and by inserting the last integral 
representation (4.1) for M 2 (a). We thus obtain 

r.TT/2 



y rdrf(r) = -J^ a 2 M 2 {a) da = - J cos 3 d0 = - , 

/•oo d /"t/2 Z" 71 "/ 2 /-oo t-> 

/ rdrS 2 (r) = - cosf3df3 cos0d0/ e (i?cos /?-«;)/ cos ^ = 

J.R ^ JO io J R cos 13 4 



while the second contribution to Si(r) is negligible. 
These exact expressions lead to the estimate 

ln# eff ^-A (4.10) 

It turns out that this leading-order estimate has an interesting correction of relative 
order R\nR, originating in the space integral of the last term in the left-hand side of 
eq. (4.6), which can be estimated as follows. The R — > limit F (r) of the function F(r) 
obeys the equation 

F ° (r) "/n Fo(/9) = /(r) "3^ : / M 2 (a)da. (4.11) 
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We observe that the right-hand side of this equation diverges as — 4/(37rr) as r — > 0. So 
does its solution: 

F (r)^-^- (r-0), (4.12) 

since the integral in the left-hand side of eq. (4.11) is less singular. The space integral of 
the last term in the left-hand side of eq. (4.6) can be estimated by setting 5 = a — s ps 
a + Rcos/3 - r > 0: 

00 //■ r , / , AR [°° ^ f°° ^ . . n r /2 cos/3d/3 



jf rdr JJf ( P )^-^ J drj^ M 2 (a) da J 



(R 2 + S 2 -2RS cos /3) 



1/2 



4# f 00 „ , , , / 2(7 7T \ 4i? / R TV 

M 2 {a) da In — + - + 1 « — In - + 1e - 



3n 2 J zv ' \ R 2 J 3tt 2 V 4 2 

(4.13) 

with 7 B being Euler's constant. 

We thus obtain the more accurate estimate 

i Reft 4£ 16 , R , 

In— — « H -In— + lna, (4.14) 

£ 3# 3tt 2 £ ' v ; 

i.e., 

^eff^a£^ 37r2 e"3i?, (4.15) 

in physical units. It is worth noticing that the correction in RlnR derived in eq. (4.13) 
is responsible for the occurrence of the power-law prefactor in the result (4.15). The 
dimensionless absolute prefactor a receives contributions from several corrections of relative 
order R to the above estimates. A fit of the numerical data shown in Figure 5 leads to 
a « 1.20. 

The capacitance Q thus reads, from eq. (2.7), 

_ 3R ( 3 . R 9 . r \ R 2 f . 

Q~ h[-oln In— — . (4.16) 

^ U \ix 2 i 16 at) £ 2 K ' 

The leading behavior of the capacitance is therefore an intrinsic characteristic of the ab- 
sorbing cylinder, while the length scale r , involving the geometry of the sample, only 
appears logarithmically in the correction term. 

• Large-radius behavior 

If the radius of the cylinder is large (R^$> 1), eq. (4.3) again becomes the SM equa- 
tion (3.22) of RTT in a half-space with a free boundary. By requiring that the asymptotic 
behavior (3.24) of its solution is proportional to the solution of the diffusion equation, 
J(r) = ln(r/-R e ff), in the range 1 C r < J? [see eq. (4.5)], we obtain the estimate 

R eS » Re~^/ R ^Rfl-^\, (4.17) 
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in physical units, for R ^> £, where To is again Milne's extrapolation length (3.25). This 
result agrees, to leading order, with the prediction (2.21) of the diffusion approximation, 

Figure 5 shows a logarithmic plot of the size factor R eG /R*g against R/L The 
numerical data interpolate between the asymptotic behaviors (4.15) and (4.17), listed in 
Table 2. 

4.3 The polarizability of a transparent cylinder 

In analogy with eq. (3.28), the RTT equation for a transparent cylinder reads 



J(x) = JJl(x-na T ) + JJ J(x-ikt), (4.18) 

with 

a T = a + 2Rcos/3 = a + 2(R 2 - r 2 sin 2 a) 1/2 . (4.19) 

In order to calculate the polarizability P, along the lines of sections 2.3.2 and 3.3, we 
look for a solution to eq. (4.18) of the form 

J(x) = xK(r) = rK(r) cos 9, (4.20) 

with K{r) going to unity at large distances. Eq. (4.18) becomes 

rK(r)= // (r - a T cos a)K(p T ) + // (r - acosa)K(p), (4.21) 



where p has been defined in eq. (4.4), and with 

p T = ( r 2 + of - 2ra T cosa) 1/2 . (4.22) 

We set 

K(r) = l + ^ + G{r), (4.23) 

where G(r) again represents the short-ranged skin-layer correction to the solution of the 
diffusion equation. The equation for G(r) reads 



rG(r)+ // (a cos a-r)G(p)+ // [(a T cos a-r)G{p T )+{r-a cos a)G(p)] =S 1 (r)+PS 2 (r), 

i/full JJl 

(4.24) 

with 

^i(^) = // ip — &t) cos a = — 2P // cos a cos /3, 

/ x (4-25) 

^ . . . . // it — <7t cos a a cos a — r x 

«,(r) = -,(r) + jf(— ^ + — -3- 
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and 

/ s 1 ff r — a cos a 1 f°° , , . . , , , 

g(r) = -- - 2 = - / M 2 (a) da. (4.26) 

r i/full P r 

• Small-radius behavior 

The behavior of the polarizability P for J? C 1 can again be derived by taking the 
Fourier transform at q — > of eq. (4.24), after multiplication by cos6>. The Fourier 
transform of the function J(x) defined in eq. (4.20) reads 

X(q) = 2z7r/ r 2 dr K^J^qr) ^ inq r s drK(r) (q -> 0), (4.27) 
io io 

where Ji is the Bessel function. Therefore both sides of eq. (4.24) have to be integrated 
with a weight r 2 dr. The integral of S\(r) can again be estimated by changing variables 
from a to P and from r to w. We thus obtain 

/ r 2 drS 1 (r) = — — / cos 2 [3 d[3 cos0d<W e {R <*»(>-*>), >co S <j> w dw 

JR 71 Jo Jo J R cos (3 

R 2 



= - — (l + R), 



•7T/2 2 

3' 



/ r 2 drg(r) = - a 2 M 2 (a) da = / cos 3 </>d</> 
io ^ io io 

(4.28) 

while the second contribution to £2 (r) is negligible. This leads to the estimate 

P«-^fl + ^V (4.29) 



2 V 

in physical units, up to a correction of relative order R 2 /£ 2 . 

• Large-radius behavior 

When the radius of the sphere is large 1), eq. (4.21) again becomes a RTT 

equation in a half-space. With the definitions (3.20), (3.38), and setting again K{r) = T(r), 
we obtain 

r(r)- / dr , [M 1 (r-r , ) + M 1 L (r + r , )]r(r / ) = 0, (4.30) 
Jo 

where 

M 1 (x)= [[ 5(x-acosa)= f ^e" 1 ^, 

iifull i() 2 P ^ g-j^-j 

M[(x)= ff (l-2cos 2 a)5(x-acosa) = [ ^(1 - 2 / u) e - |x|/ '\ 

i/full io 2 P 
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We thus obtain the SM equation in a half-space with partial reflection at the boundary, 
characterized by a reflection coefficient 

R(fj) = 1 - 2fi, (4.32) 

which is again negative for a small enough, namely [i = cos a > 1/2, i.e., a < tt/3. 
The normalized solution to eq. (4.30) again behaves as 

T(t)^t + t 01 (t>1). (4.33) 

We have evaluated tqi numerically, again along the lines of section V of ref. [12], obtaining 
thus 

r i « 0.2137. (4.34) 

By requiring that expression (4.33) is proportional to the solution of the diffusion 
equation, K(r) = 1 + P/r 2 , in the range 1 C r <C i? [see eq. (4.23)], we obtain the 
estimate 

R 3 _. 02 ^ 2r i£ 



. r -R z 1 - =^ , (4.35) 

i2 + 2r i£ V # / 

in physical units, for R^> £. This result agrees, to leading order, with the prediction (2.28) 
of the diffusion approximation, P dlff = — R 2 . 

Figure 6 shows a plot of the size factor _p/_p dlff against R/£. The numerical data 
interpolate between the asymptotic behaviors (4.29) and (4.35), listed in Table 2. 

4.4 The polarizability of a reflecting cylinder 

In analogy with eq. (3.46), the RTT equation for a reflecting cylinder reads 

J(x) = JJAyr) + jj A* ~ n<r), (4.36) 

with 

y a (x, n,s):\ XR = r ~ S ° C ° S ? + ^ ~ f )™^ a) ' (4.37) 

JHK ' ' ' \ y R = s sin a + (a - s )sm(2/3 - a). v ' 

We again look for a solution to eq. (4.36) of the form 

J(x) = xK(r) = rK(r) cos 0, (4.38) 
with K(r) going to unity at large distances. Eq. (4.36) becomes 

rK(r) = JJ^ x R K(p R ) + JJ (r - a cos a)K{p), (4.39) 

with the notation (3.49). Setting 

K{r) = l + ^ + G{r), (4.40) 
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the equation for G(r) reads 
rG(r)+ !! (acosa-r)G(p)+ [[ [(r-acos a)G(p)-x R G(p R )] = S 1 (r)+PS 2 (r), (4.41) 

i/full JJl 

with 



Si(r) = 2 11 (o" — So) cos j3 cos(ct — /3), 
ff f x R a cos g - r 

where y(r) has been defined in eq. (4.26). 

• Small-radius behavior 

The behavior of the polarizability P for R <C 1 can again be derived along the lines 
of the previous sections. We have 

r 2 dr5'i(r) = — / cos 2 (3d(3 / cos 2 </>d^ 

R 77 JO JO 

poo 

x / e (i?cos/3 - w)/cos ^(«;cos/3 + Psin 2 /3)d«; (4.43) 

J R cos /3 



-(Mr 

hence the estimate 

^ 3P£ 

^~-g-, (4.44) 
in physical units, up to a correction of relative order Rjt. 

• Large-radius behavior 

In analogy with section 3.5, we have been led to study the solution to the RTT 
equation, including corrections of order 1/R 3 included. Setting 

we obtain successively 

r 1 (r) = -r, r 2 (r) = ^r 2 , T 3 (r) = -2r 3 - jjr + 73 (r). (4.46) 

By requiring that expression (4.45) is proportional to the solution of the diffusion 
equation, K(r) = 1 + P/r 2 , in the range 1 « r « E [see eq. (4.40)], we obtain the 
estimate 

/ 4/2 \ 

P " R2 1+ 5JP ' (447) 
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in physical units, for R 3> £. This result agrees, to leading order, with the prediction (2.29) 
of the diffusion approximation, P dlff = R 2 . 

Figure 7 shows a plot of the size factor _p/_p dlff against R/t. The numerical data 
interpolate between the asymptotic behaviors (4.44) and (4.47), listed in Table 2. 

We have also checked that the alternative approach sketched at the end of section 3.5 
permits to recover the result (4.47). 

5 DISCUSSION 

In this paper, we have first recalled that the long-range effects of an isolated embed- 
ded object on the diffusive propagation of radiation in an otherwise homogeneous, non- 
absorbing turbid medium can be described within the diffusion approximation [2, 4-7]. 
The leading term of a multipole expansion involve a (negative) charge q (for an absorbing 
object) or a dipole moment p (for a non-absorbing object). We have emphasized the role 
of the associated intrinsic characteristics of a spherical object, namely its capacitance Q 
and its polarizability P. For cylinders, we have underlined that the right intrinsic quan- 
tity to be considered is its effective radius -R e ff- The intrinsic quantities Q, R e s, P are 
related to the observables q and p, which depend on the geometry of the sample, through 
eqs. (2.4), (2.7). We have shown that RTT provides the adequate framework to evaluate 
the characteristics Q, _R e ff, and P of mesoscopic objects, with a radius R comparable to 
the mean free path £, the latter being much larger than the wavelength A. Considering for 
simplicity the case of isotropic scattering of scalar waves, we have worked out this program 
for spheres and cylinders of radius R, in several situations: 

(i) totally absorbing (black) object, 

(ii) transparent object, with the same optical index as the medium, 

(iii) totally reflecting object. 

For large objects (R 3> t), we recover the predictions of the diffusion approxima- 
tion [2, 6], and we give the analytical form of the first correction to this approximation, 
which is either proportional to tjR or to (£/R) 2 . In the opposite regime of small objects 
(R <C £), our results have some overlap with those obtained by the partial- wave expansion 
of Lancaster and Nieuwenhuizen [7]. In the intermediate situation, i.e., for generic values 
of the ratio R/£, the characteristics Q, R e g, or P differ from the predictions Q dlfT , R^g, 
or _P dlff of the diffusion approximation by a multiplicative size factor, which depends con- 
tinuously on R/£, and smoothly interpolates between the analytical results pertaining to 
both limiting regimes. For cases (i) and (ii), the size factor is a monotonically increas- 
ing function of R/£, converging to unity from below, with a correction of order £/R. For 
case (iii), the size factor is a decreasing function of R/i, converging to unity from above, 
with a correction of order (£/R) 2 . 

The characteristics Q, R e s, and P of spheres and cylinders are expected not to be fully 
universal, but rather to weakly depend on microscopic details of the scattering process. 
The associated size factors have been calculated in this work for isotropic scattering of 
scalar waves. They would be slightly different in other situations, such as anisotropic 
scattering of scalar waves, or Rayleigh scattering of vector waves. This weak dependence 
can be illustrated on the example of large absorbing objects. Indeed, for both spheres 
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and cylinders, the first correction to the diffusion approximation has been shown to be 
proportional to the extrapolation length To of the problem of RTT in a half-space with a 
free boundary. This quantity is known to be tq ~ 0.718211 for very anisotropic scattering 
of scalar waves [13], in units of the scattering mean free path £*, and tq ~ 0.712110 
for Rayleigh scattering of electromagnetic waves [14], these numbers being respectively 
1.09% and 0.23% above the Milne value (3.25), To m 0.710446, corresponding to isotropic 
scattering of scalar waves. 

Finally, in all the situations considered in this work, the effect of the embedded object 
on a light ray can be described in simple terms. In cases (i) and (iii) radiation does not 
enter the object at all, while in case (ii) it propagates ballistically inside the object. The 
present RTT approach can be generalized, at least in principle, to mesoscopic spheres or 
cylinders with yet other kinds of optical properties. For instance, situations where radiation 
is partly absorbed in the embedded object and/or in the medium, can be dealt with rather 
simply, in the general case where the corresponding absorption lengths are comparable 
to the size of the object. A more complex situation is e.g. that of a transparent object, 
with an optical index different from that of the medium, so that light undergoes multiple 
reflections inside the object. In this case the RTT equation receives an infinite series of 
contributions from the multiply reflected light rays, weighted by the appropriate extinction 
factors. Furthermore, as opposed to the cases treated in the present paper, the light rays 
are not contained in a plane in general, so that the relevant integral equations keep two 
angular variables. 

Acknowledgments 

Th.M.N. acknowledges useful discussions with G. Maret, R. Maynard, M. Nieto- 
Vesperinas, and J. Ricka. 

This work has been supported by the HCM network entitled Transport of Light in 
Strongly Scattering Media, under contract no. ERBCHRXCT 930373. 



29 



REFERENCES 



[I] M.C.W. van Rossum and Th.M. Nieuwenhuizen, Rev. Mod. Phys. (to appear in 
January 1999). 

[2] P.N. den Outer, Th.M. Nieuwenhuizen, and A. Lagendijk, J. Opt. Soc. Am. A 10, 
1209 (1993). 

[3] J.C. Schotland, J.C. Haselgrove, and J.S. Leigh, Appl. Opt. 32, 448 (1993). 

[4] S. Feng, F.A. Zeng, and B. Chance, Appl. Opt. 34, 3826 (1995). 

[5] X.D. Zhu, S.P. Wei, S. Feng, and B. Chance, J. Opt. Soc. Am. A 13, 494 (1996). 

[6] J.C.J. Paasschens and G.W. 't Hooft, J. Opt. Soc. Am. A 15, 1797 (1998); J.C.J. 
Paasschens, Ph.D. thesis (Leiden, 1997). 

[7] D. Lancaster and Th.M. Nieuwenhuizen, Physica A 256, 417 (1998). 
[8] S. Chandrasekhar, Radiative Transfer (Dover, New- York, 1960). 

[9] V.V. Sobolev, A Treatise on Radiative Transfer (Van Nostrand, Princeton, N.J., 1963). 

[10] A. Ishimaru, Wave Propagation and Scattering in Random Media, in 2 volumes (Aca- 
demic, New- York, 1978). 

[II] H.C. van de Hulst, Multiple Light Scattering, in 2 volumes (Academic, New- York, 
1980). 

[12] Th.M. Nieuwenhuizen and J.M. Luck, Phys. Rev. E 48, 569 (1993). 
[13] E. Amic, J.M. Luck, and Th.M. Nieuwenhuizen, J. Phys. A 29, 4915 (1996). 
[14] E. Amic, J.M. Luck, and Th.M. Nieuwenhuizen, J. Phys. I France 7, 445 (1997). 
[15] J.C.J. Paasschens, private discussion. 



30 



CAPTIONS OF FIGURES AND TABLES 



Figure 1: Geometrical characteristics of rays which hit the object under study (a sphere 
centered at O in section 3, or a cylinder whose axis is parallel to the z-axis and passes 
through O in section 4) . Table 1 contains the expressions of various quantities attached to 
this Figure. 

Figure 2: Plot of the size factor Q/Q dlS of the capacitance of an absorbing sphere, against 
the ratio R/t. Full line: outcome of the numerical analysis explained in the text. Dashed 
lines: small-radius and large-radius behaviors, as listed in Table 2. 

Figure 3: Same as Figure 2, for the size factor p/p dlff of the polarizability of a transparent 
sphere. 

Figure 4: Same as Figure 2, for the size factor _p/_p dlff f the polarizability of a reflecting 
sphere. 

Figure 5: Same as Figure 2, for the size factor R e g/R^ of the effective radius of an 
absorbing cylinder (logarithmic scale). 

Figure 6: Same as Figure 2, for the size factor p/p dlff of the polarizability of a transparent 
cylinder. 

Figure 7: Same as Figure 2, for the size factor _p/_p dlff of the polarizability of a reflecting 
cylinder. 

Table 1: Summary of useful definitions, including distances and other quantities attached 
to the rays constructed in Figure 1, domains I and II in both spherical and cylindrical 
geometries, and the associated integration measures. 

Table 2: Asymptotic behavior at small and large radius of the various characteristics 
of embedded spherical and cylindrical objects studied in the text. Absorbing objects are 
characterized by their capacitance Q or their effective radius R e ff, while non-absorbing, 
either transparent or reflecting, objects are characterized by their polarizability P. The 
large-radius behaviors agree, to leading order, with the predictions of the diffusion approx- 
imation, recalled in section 2.3. 
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Table 1 



Geometric quantities 



r = \OX\ 

w = \AX\ = rcosa = (r 2 - R 2 sin 2 (3) 1/2 

\OA\ = rsina = Rsin/3 

\AB\ = \AC\ =Rcos(3= (R 2 - r 2 sin 2 a) 1 / 2 

sinao = R/r 

s = \XC\ = w — Rcos f3 

r>R, so > 



Spherical geometry 
(Section 3) 



I : < a < ao and so < s < oo 

< a < ao and < s < so 

II : < or 

ao < a < it and < s < oo 



// - = / 



sin a da 



-e" s ds--- 



cos a 



Cylindrical geometry 
(Section 4) 



I : < a < ao and so < o < oo 

{0 < a < ao and < a < so 
or 
ao < a < tt and < a < oo 

// ...= // — M 2 (o-) do-' • • 

i//,//,full i//,//,full 7r 

1 M 2 (a) d<7 



// 



i2 



7rr 



r/2 



cos /3 d/3 



./o 



cos a 
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Table 2 



Object 


Small-radius behavior 

(R<^£) 


Large-radius behavior 

(R > I) 


absorbing sphere 
(Sec. 3.2, Fig. 2) 






m 2 
u 




r « 0.710446 


transparent sphere 
(Sec. 3.3, Fig. 3) 




3 \ 




) 

/ 


pdiff 

r i ~ 0.4675 


reflecting sphere 
(Sec. 3.4, Fig. 4) 




P ft; 


4 




„ P 3 / 6£ 2 \ 

pdiff 


absorbing cylinder 
(Sec. 4.2, Fig. 5) 

V 7 0/ 


R e fl ft 


-(? 


16 

J 6 


4£ 
3~R 


P eff ^ ^ (l " ^) 

pdiff 
^eff 






a ss 1.20 




r « 0.710446 


transparent cylinder 
(Sec. 4.3, Fig. 6) 


P ( 


P 2 



2 




1 


pdiff 

T i « 0.2137 


reflecting cylinder 
(Sec. 4.4, Fig. 7) 




P ftj 


mi 

8 




pdiff 
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